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ABSTRACT 

The  sound  generated  by  rainfall  at  sea  is  caused  by  raindrops  of  a  wide  range  of  sizes 
and  angles  of  incidence,  which  fall  at  their  terminal  velocities.  The  purpose  of  this  labo- 
ratory research  has  been  to  make  complete  acoustical  measurements  of  the  sound  gen- 
erated by  single  water  drops  striking  the  water  surface  at  their  terminal  velocities  for 
normal  and  oblique  incidence.  These  measurements  have  included  the  total  acoustic 
energy,  peak  axial  pressure,  frequency  spectrum  and  polar  radiation  pattern. 

Depending  on  the  drop  size  and  the  angle  of  incidence,  many  drops  falling  at  their 
terminal  velocities  create  bubbles.  At  all  angles  of  incidence  studied  here  the  sound  ra- 
diated by  an  individual  bubble  contains  more  energy  than  the  sound  from  an  individual 
impact. 

These  results,  using  terminal  velocities  and  oblique  trajectories,  are  very  different 
from  the  published  normal  incidence,  non-terminal  velocity  characterizations.  For  ex- 
ample bubble  frequencies  other  than  the  well  known  14  kHz  peak  are  found.  Also  the 
energy  of  the  impact  sound  increases  significantly  for  larger  drop  sizes  and  for  larger 
deviations  from  normal  trajectories.  Furthermore,  drops  of  diameter  0.8  mm  to  1  mm, 
which  always  produce  bubbles  at  normal  incidence,  create  bubbles  only  about  10  %  of 
the  impacts  at  oblique  incidence.  These  observations  provide  specific  reasons  for  the 
previously  unexplained  broadening,  shifting  and  reduction  in  magnitude  of  the  14  kHz 
spectral  peak  of  the  rain  noise  in  the  presence  of  winds  at  sea. 


in 


/JJ&0 


e/i 


TABLE  OF  CONTENTS 

I.  INTRODUCTION    

II.  THEORY  OF  THE  EXPERIMENT    5 

III.  EXPERIMENT   

A.     EQUIPMENT     

7 

1.  Mechanical   

7 

2.  Optical    

7 

3.  Acoustical     

8 

4.  Electronic    

5.  Experimental  Procedure     

B.  NORMAL  INCIDENCE     10 

.  12 

1      Peak  Axial  Pressure     

20 

2.  Total  Acoustic  Energy     

26 

3.  Frequency  Specta    

...  30 

4.  Radiation  Pattern    

C.  OBLIQUE  INCIDENCE     34 

"J  A 

1.  Peak  Axial  Pressure     

"3  -j 

2.  Total  Acoustic  Energy     

46 

3.  Frequency  Spectra     

A      Percentacc  of  Drops  Creating  Bubbles   

,.    .  "  n  56 

5.     Radiation  Pattern    

59 
IV.     CONCLUSIONS 

APPENDIX  A.     PREAMPLIFIER  1     6l 

APPENDIX  B.     PREAMPLIFIER  2    62 

LIST  OF  REFERENCES    6j 


IV 


INITIAL  DISTRIBUTION  LIST 


64 


LIST  OF  FIGURES 

Figure     1.  Noise  at  sea     2 

Figure     2.  Regions  of  regular  and  random  bubble  creation    3 

Figure     3.  Experimental  setup    9 

Figure     4.  Bubble  creation  region  including  small  drops    11 

Figure     5.  Acoustic  pressure  of  an  impact  and  a  bubble 13 

Figure     6.  Acoustic  pressure  caused  by  impact    14 

Figure     7.  Acoustic  pressure  caused  by  bubble    15 

Figure     S.  Average  peak  axial  pressures  for  impacts    16 

Figure     9.  Average  peak  axial  pressures  for  bubbles    17 

Figure   10.  Peak  axial  pressure  histograms  of  impacts     18 

Figure   11.  Peak  axial  pressure  histograms  of  bubbles     19 

Figure  12.  Average  acoustic  energies  radiated  by  the  impacts    22 

Figure  13.  Average  acoustic  energies  radiated  by  the  bubbles     23 

Figure   14.  Acoustic  energy  histograms  of  the  impacts   24 

Figure  15.  Acoustic  energy  histograms  of  the  bubbles 25 

Figure   16.  Frequency  spectrum  of  the  impact  sound 27 

Figure   17.  Peak  frequencies  of  impacts   28 

Figure   IS.  Resonance  frequencies  of  bubbles     29 

Figure  19.  Fight  hydrophone  array  with  20  degere  separations     31 

Figure  20.  Eight  hydrophone  array  with  10  degree  separations     32 

Figure  21.  Polar  radiation  pattern  of  normal  incidence  drops     33 

Figure  22.  Average  peak  axial  pressures  for  impacts    35 

Figure  23.  Average  peak  axial  pressures  for  bubbles    36 

Figure  24.  Max.  min.  and  avg.  impact  acoustic  energies  for  0.83  mm  drops    38 

Figure  25.  Max.  min.  and  avg.  impact  acoustic  energies  for  0.914  mm  drops    39 

Figure  26.  Max.  min.  and  avg.  impact  acoustic  energies  for  0.9S5  mm  drops    40 

Figure  27.  Average  impact  acoustic  energies  of  the  drops    41 

Figure  28.  Radiated  bubble  energy  for  0.83  mm  drops     42 

Figure  29.  Radiated  bubble  energy  for  0.914  mm  drops     43 

Figure  30.  Radiated  bubble  energy  for  0.985  mm  drops     44 

Figure  31.  Average  bubble  acoustic  energies  for  oblique  incidence     45 


VI 


Figure  32.  Frequency  spectrum  of  an  impact  sound     47 

Figure  33.  Avg.  peak  frequencies  of  impact  sound  for  oblique  inc.  drops    48 

Figure  34.  Max.  min.  and  avg.  resonance  frequencies  of  bubbles  for  0.83  mm    ....  49 

Figure  35.  Max.  min.  and  avg.  resonance  frequencies  of  bubbles  for  0.914  mm     ...  50 

Figure  36.  Max.  min.  and  avg.  resonance  frequencies  of  bubbles  for  0.9S5  mm     ...  51 

Figure  37.  Avg.  resonance  frequencies  of  bubbles  for  oblique  incidence  drops    ....  52 

Figure  38.  Percentage  of  drops  that  create  bubbles    54 

Figure  39.  The  ratio  of  the  average  bubble  energy  to  the  impact  energy     55 

Figure  40.  Polar  radiation  pattern  of  bubbles  for  10  degrees  incidence 57 

Figure  41.  Polar  radiation  pattern  of  bubbles  for  20  degrees  incidence 58 


vu 


LIST  OF  SYMBOLS  AND  ABBREVIATIONS 


SYMBOL 

UMTS 

NAME/  DESCRIPTION 

P 

kg,  m3 

Water  density 

c 

m/  sec 

Speed  of  sound  in  water 

w 

ml  sec 

Particle  velocity 

R 

meter 

Range  of  hydrophone  on  the  axis 

P 

pascals 

Pressure 

k 

radians' 

meters 

Acoustic  wave  number 

KE 

Joules 

Kinetic  energy 

m 

kg 

Mass 

e 

degrees 

Angle  from  the  vertical  acoustic  axis 

V 

meters,' 

sec 

Velocitv 

vm 


ACKNOWLEDGEMENTS 

I  wish  to  express  my  sincere  appreciation  to  Professor  Herman  Medwin  of  the 
Physics  department,  U.S  Naval  Postgraduate  School,  for  his  guidance,  encouragement 
and  patience  which  made  this  study  a  great  learning  experience;  to  Professor  Jeffrey  A. 
Nystuen  for  his  invaluable  help  in  all  parts  of  the  experiment,  to  LTJG  Olmez  for 
spending  his  summer  vacation  in  my  laboratory,  trying  hard  for  video  film  of  the  drops 
and  to  my  wife  Fazilet  for  her  support,  understanding  and  love. 


IX 


I.     INTRODUCTION 

There  have  been  no  theoretical  studies  of  the  sound  generated  by  rainfall  at  oblique 
incidence,  which  most  probably  always  occurs  because  of  the  wind  as  well  as  the 
roughness  of  the  sea  surface.  Instead  there  have  been  theoretical  studies  (Franz,  1959), 
(Longuet-Higgins,  19S9),  (Oguz  and  Prosperetti,  1989)  for  normal  incidence.  Also  there 
have  been  no  experimental  results  published  before  for  oblique  incidence  of  single  drops. 
Instead  there  have  been  several  experiments  for  drops  at  normal  and  non-terminal  ve- 
locities. The  purpose  of  this  thesis  is  to  study  the  sound  generated  by  raindrops  striking 
the  water  surface  at  oblique  incidence  and  at  normal  incidence. 

Several  previous  experiments,  both  with  natural  and  artificial  rain,  have  verified  that 
the  main  feature  of  the  noise  spectrum  generated  by  rain  at  sea  is  a  peak  at  about  14 
kHz.  (Nystuen  and  Farmer,  1987).  This  peak  broadens,  shifts  to  higher  frequencies  and 
decreases  in  magnitude  as  the  wind  speed  increases  (Figure  1). 

G.J. Franz  (1959)  concluded  a  laboratory  study  of  single  drops  by  asserting  that  rain 
drops  create  noise  by  two  different  mechanisms;  IMPACT  and  OSCILLATION  of  the 
bubbles  created  by  the  drops.  He  also  concluded  that  at  low  impact  velocities  bubble 
noise  dominates  but  at  high  impact  velocities  bubbles  are  less  likely  to  be  formed  so  that 
the  sound  generated  by  a  raindrop  at  its  terminal  velocity  could  be  considered  as  the 
sound  from  the  impact  only. 

Later  studies  showed  that  oscillating  air  bubbles  are  more  important  noise  sources 
than  Franz  had  argued.  Pumphrey  and  Crum  (1989)  conducted  experiments  for 
normally  incident  single  drops  with  different  sizes  and  different  impact  velocities  as 
shown  in  Figure  2.  Their  explanation  follows  :  the  shaded  area  is  a  region  of  consistent 
bubble  creation.  In  the  range  below  this  region  the  impact  is  not  energetic  enough  to 
create  a  bubble.  Above  the  region  the  energy  is  excessive.  At  still  higher  velocities  and 
drop  sizes  (upper  right  corner)  there  is  another  region  for  which  bubbles  occur  in  an  ir- 
regular and  random  fashion.  The  curve  at  the  left  of  Figure  2  shows  the  terminal  ve- 
locity of  drops  as  a  function  drop  diameter.  The  range  of  drop  diameters  where  this 
terminal  velocity  curve  intersects  the  shaded  region  and  causes  bubble  entrainment  is  the 
narrow  region,  stated  by  them  as 

O.S  mm  <  DROP  DIAMETER  <  1.1  mm 


Figure  1.      Noise  at  sea:     This  figure  is  taken  from  Nystuen  and  Farmer  (1987). 
It  shows  the  noise  spectra  at  sea  for  different  wind  conditions. 


Figure  2.  Regions  of  regular  and  random  bubble  creation:  This  figure  is  taken 
from  Crum  and  Pumphrey  (1989).  The  shaded  area  shows  the  regular 
bubble  creation  area  and  the  upper  right  corner  is  the  irregular  (random) 
bubble  creation  area. 


Crum  and  Pumphrcy  argue  that  the  bubbles  entrained  by  these  drops  have  narrowly 
defined  sizes  and  that  they  oscillate  at  a  frequency  around  14.5  kHz.  Only  drops  in  this 
range  entrain  bubbles  with  high  probability.  Since  bubble  sound  is  louder  and  in  a 
narrower  frequency  band  than  the  impact  sound,  they  conclude  that  these  particular 
drops  lead  to  the  formation  of  the  14  kHz  peak  in  rainfall.  Furthermore,  since  drops 
of  this  bubble-producing  size  are  likely  present  in  both  natural  and  artificial  rain  they 
conclude  that,  even  though  a  small  fraction  of  all  drops  are  capable  of  creating  bubbles, 
the  acoustical  effect  of  this  range  of  drops  is  dominant. 

In  this  thesis  previous  experimental  results  for  normal  incidence  are  first  verified  for 
small  drop  sizes  which  strike  the  water  surface  at  their  terminal  velocities.  Then  acous- 
tical radiation  from  single  water  drops  striking  the  water  surface  at  terminal  velocity  and 
oblique  incidence  is  studied. 


II.     THEORY  OF  THE  EXPERIMENT 

Since  the  radiation  pattern  of  the  noise  from  the  bubble  as  well  as  the  impact  noise 
is  assumed  to  be  dipole  radiation  for  both  normal  and  oblique  incidence  we  can  find  the 
total  acoustic  energy  produced  by  knowing  the  axial  pressure.  Assuming  pex(r)  is  the 
measured  acoustic  pressure  on  the  axis  of  the  dipole  at  range  R,  px  is  the  pressure  at  1 
meter,  pUx  is  the  acoustic  pressure  on  the  axis  at  1  meter.  For  a  dipole  source 
P\  =  P\j,  cos  $  i  where  6  is  the  angle  from  the  axis.  So  the  acoustic  energy  can  be  written 
as  : 


Energy  = 


{INTESSITY)dSdt 


Energy  = 


p^dSdt 


rs 


where  vv,  is  the  particle  velocity  and  dS  is  the  element  of  area. 

P\ 

When  kR  >   3  we  can  use  w,  =  —tt  with  less  than  10  %  error  in  magnitude,  where 

p  is  the  density  of  the  water,  c  is  the  speed  of  sound  in  the  water  and  k  is  the  acoustic 
wave  number.  Doing  the  area  integration  by  using  the  ring  element  (2nR  sin  6)Rd6  we 
can  find  that 


Energy 
Time 


{In 


EL 

pc 


'f 


7  ' 

plax  cos'tf  sin  Odd 


=  (2  —  )RPlax 

Because  the  impact  sound  and  the  bubble  sound  are  not  continuous,  it  is  convenient  to 
do  the  integration  over  time  bv  a  summation  for  elements  of  time  in  the  digitized  form. 


Energy  =  {In -j^)  J  p]axt 


This  can  be  compared  with  the  kinetic  energy  of  the  rain  drop  just  before  the  impact 
time  where 

KE  =  -r  mv 


III.     EXPERIMENT 

A.     EQUIPMENT 

1.  Mechanical 

An  Eppendorf  digital  pipette  which  can  produce  accurately  metered  drops  from 
0.1  n  liter  to  10  \x  liter  (0.58  mm  to  2.67  mm  diameter)  was  used  to  create  definite  sizes 
of  artificial  raindrops.  The  dropper  was  placed  on  a  plate  located  at  a  fixed  position  of 
1.8  meter  over  the  water  surface  to  make  sure  that  all  drops  strike  at  the  same  point  over 
the  axial  hydrophone.  At  this  height  all  drops  smaller  than  0.985  mm  diameter  reach 
their  terminal  velocities  (Wang  and  Pruppacher,  1977).  The  hydrophones  were  placed 
in  an  anechoic  water  tank  10  ft  x  10  ft  x  9.5  ft  deep  with  redwood  wedge  absorbers. 

Oblique  incidence  was  obtained  by  using  a  fan  to  blow  the  drop  horizontally, 
immediate  after  it  started  to  fall.  The  water  surface  was  not  disturbed  by  the  wind. 

2.  Optical 

To  observe  the  angle  of  incidence  for  the  oblique  incidence  case  a  video  camera 
with  7  lux  light  requirement  and  high  speed  shutter  was  used  together  with  an  HQ  video 
player.  The  player  (Hitachi  VT  498  EM)  had  adjustible  time,  frame  by  frame  slow  mo- 
tion and  view  freezing  capability  which  made  possible  both  drop  velocity  measurements 
as  well  as  incidence  angle. 

Alternatively  a  camera  (Canon  AE-1)  was  also  used  together  with  an  adjustable 
frequency  stroboscope  to  observe  both  the  impact  velocity  and  the  incidence  angle  of  the 
drops.  For  both  techniques  we  used  a  5  cm  x  5  cm  string  grid  behind  the  paths  of  the 
drops  to  locate  the  position  of  the  drop. 

3.  Acoustical 

An  ITC  6050  C  hydrophone  with  high  sensitivity  (-163  dB  ref  IV  per  /j.  Pa)  and 
flat  frequency  response  from  1  to  40  kHz  was  used  to  observe  the  impact  and  bubble 
noise  for  the  normal  incidence  case.  For  oblique  incidence  a  nine  hydrophone  array  was 
used.  These  hydrophones  were  either  LC10  of  sensitivity  -198  dB  or  our  own 
hydrophones  constructed  of  6  mm  diameter  cylindrical  lead  zirconatc  elements.  The 
sensitivities  of  our  own  hydrophones  are  -205  dB.  The  array  was  mounted  on  a  J  8  inch 
diameter  thin  semi-circular  rod  bent  to  30  centimeter  radius.  The  hydrophones  were 


placed  even'  10  or  20  degrees  along  the  rod.  The  hydrophones  were  calibrated  by  ab- 
solute reciprocity  and  by  inter-comparison,  using  a  ping  or  the  impulse  from  a  reference 
raindrop. 

4.  Electronic 

To  investigate  the  polar  radiation  pattern,  nine  hydrophones,  nine  preamplifiers, 
nine  high  pass  filters,  were  the  input  to  an  IBM  PC  -XT  computer  with  an 
R"C-Electronics  analog  to  digital  converter  board.  A  software  package,  called 
COMPLTERSCOPE,  having  a  maximum  AD  conversion  rate  of  106  samples  per  sec- 
ond and  12  bit  amplitude  resolution  processed  the  signals.  We  made  our  own  pream- 
plifiers as  shown  in  Appendices  1  and  2. 

Also  to  calibrate  the  hydrophones,  an  HP-4192  Impedance  Analyzer,  a  signal 
generator,  a  power  amplifier  and  a  transducer  were  used.  The  transducer  was  placed  in 
the  far  field  and  was  driven  with  pulsed  signals  to  prevent  scattering  from  the  tank  walls. 

5.  Experimental  Procedure 

The  experimental  setup  was  as  shown  in  Figure  3.  The  noise  created  by  the 
bubble  and  the  impact  was  recorded  by  the  hydrophone  at  depth  of  1  m,  directly  below 
the  point  of  impact  for  normal  incidence  and  at  depth  of  20  cm  for  oblique  incidence  and 
for  radiation  pattern  observations.  The  output  of  the  hydrophone  was  amplified  and 
passed  through  a  high  pass  filter.  After  filtering  all  the  noise  below  8  kHz  the  signal  was 
sent  to  an  IBM  PC-XT  Computer  which  had  an  RC  Electronics  analog  to  digital  con- 
verter board  with  sixteen  input  channels. This  system  made  it  possible  to  observe  the 
outputs  of  several  different  hydrophones,  trigger  with  a  desired  voltage  level,  find  out  the 
frequency  spectrum,  maximum  pressures,  total  energies  and  radiation  patterns  of  both 
the  bubble  and  the  impact  when  used  with  the  software  package  called  Computerscope. 
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Figure  3.      Experimental  setup 


B.     NORMAL  INCIDENCE 

Initially,  different  sizes  of  drops  between  0.3  n  1  (0.83mm)  to  0.5  /j.  1  (0.985  mm  di- 
ameter) that  fell  normally  incident  to  the  tank  were  used.  The  fraction  of  drops  that 
created  bubbles  was  noted  for  all  different  drop  sizes  to  compare  with  the  experimental 
results  of  Crum  and  Pumphrey  (1989)  which  is  shown  in  Figure  2.  The  bubble  creation 
region  curve  was  re-plotted  including  our  data  points  for  near  terminal  velocities  as 
shown  in  Figure  4.  Our  data  agreed  with  theirs,  i.ev98%  of  the  drops  created  bubbles. 
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Figure  4.  Bubble  creation  region  including  small  drops:  This  figure  is  a  rcplot  of 
Figure  2  including  our  data  points,  indicated  by  symbol  x,  at  the  termi- 
nal velocitv  curve. 
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1.     Peak  Axial  Pressure 

The  noise  generated  by  a  typical  impact  is  shown  in  Figure  5  with  the  bubble 
noise  created  by  the  same  drop  17.7  ms  after  the  impact.  The  same  impact  and  bubble 
radiation  are  shown  in  Figures  6  and  7,  with  different  scales  to  see  the  details.  As  we 
can  see  from  Figure  5,  the  peak  axial  pressure  of  the  bubble  is  much  greater  than  the 
peak  axial  pressure  of  the  impact.  For  the  impact  in  Figure  6  the  peak  axial  pressure  is 
460  milivolt  after  40  dB  amplification  The  sensitivity  of  the  hydrophone  is  -163  dB  ref 
IV  per  n  pascal  and  it  is  located  20  cm  under  the  water  surface.  So  the  calculation  of 
the  peak  axial  pressure  is  as  follows  : 

-163  =  201og10^^ 


^ref 


Since  Xref  is  1  V  per  fj.  pa.  we  can  write 


ref  1 V  per  \x  Pa 


V  per  fj.  Pa 


-163       .         v 
^— -  =  log10A 


Z=7.08x  10~9 


1  6  V-Pa 

-L=  141  x  106 


V 

Pa 


=  141 


volt 


The  voltage  output  is  460  mV  after  100  times  amplification  for  the  hydrophone  is  on  the 
axis  and  20  cm  away  from  the  point  of  impact.  So  the  voltage  at  the  hydrophone  for  20 
cm  is  4.6  mV  and  for  1  m  is  920  /n  V.  To  find  the  axial  pressure  at  1  meter  we  must  simply 
multiply  this  value  by  — .  The  result  for  this  particular  case  is  0.13  pascal.  The  calcu- 
lation  of  the  peak  axial  pressure  for  the  bubbles  is  the  same  as  for  the  impacts. 

The  peak  axial  pressure  for  the  bubbles  and  the  impacts  of  different  sizes  of 
drops  were  measured. The  average  values  with  standard  deviations,  and  maximum  and 
minimum  pressures  for  different  sizes  of  drops  are  shown  in  figures  8  and  9.  The 
histograms  of  peak  axial  pressures  for  impacts  and  bubbles  are  shown  in  Figures  10  and 
11. 
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Figure  5.  Acoustic  pressure  of  an  impact  and  a  bubble.:  Impact  and  bubble  sound 
of  a  0.83  mm  diameter  drop.  The  drop  is  falling  at  its  terminal  velocity, 
perpendicular  to  the  surface.  The  time  between  vertical  grid  markers  is 
400  microseconds. 
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Figure  6.  Acoustic  pressure  caused  by  impact:  The  same  impact  shown  in  Figure 
5,  with  100  microseconds  between  vertical  grid  lines  and  the  amplitude 
increased  by  a  factor  of  2.5  to  see  the  details. 
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Figure  7.  Acoustic  pressure  caused  by  bubble:  The  same  bubble  of  Figure  5,  with 
200  microseconds  between  vertical  grid  members  and  the  amplitude  de- 
creased bv  a  factor  of  2  to  see  the  details. 
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Figure  8.  Average  peak  axial  pressures  for  impacts:  Maximum,  minimum,  aver- 
age peak  axial  pressures  and  standard  deviations  of  impacts  for  different 
drop  sizes  falling  at  normal  incidence. 
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Figure  9.  Average  peak  axial  pressures  for  bubbles:  Maximum,  minimum,  aver- 
age peak  axial  pressures  and  standard  deviations  of  bubbles  for  different 
drop  sizes  falling  at  normal  incidence. 


17 


B 

s 

QL 


30 


25 


20 


15 


10 


0 


□ 
M 


0.83  mm 

0.914  mm 
□ 

0.985  mm 


M 


□ 


□ 


t*i 


0.12  0.13  0.14  0.15 

0.125  0.135  0.145  0.155 


Pressure  (pa) 


Figure  10.  Teak  axial  pressure  histograms  of  impacts:  For  300  normal  incidence 
drops.  Since  there  were  100  of  each  size  the  vertical  axis  is  both  the 
percentage  and  the  total  numbers  of  peak  impact  pressures  observed. 
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Figure  11.  Peak  aAial  pressure  histograms  of  bubbles:  For  300  normal  incidence 
drops.  Since  there  were  100  of  each  size  the  vertical  axis  is  both  the 
percentage  and  the  total  numbers  of  peak  bubble  pressures  observed. 
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2.     Total  Acoustic  Energy 

Total  acoustic  energy  radiated  both  by  the  bubble  and  the  impact  were  calcu- 
lated by  using  the  formulas  described  in  the  previous  chapter.  The  calculation  of  the 
radiated  acoustic  energies  is  as  follows:  for  the  impact  we  measured  the  output  voltage 
at  every  2  fj.  sec  and  converted  this  value  to  pressure  using  the  formulation  described  in 
the  previous  section  of  this  chapter.  Then  by  using  the  energy  formula  given  in  Chapter 
2  (Theory  of  the  experiment),  we  found  the  acoustic  energy  radiated  by  the  impact.  We 
calculated  the  acoustic  energy  radiated  by  the  bubbles  first  again  using  the  output  volt- 
ages every  2  fx  second.  Then  for  a  simplified  calculation  we  assumed  a  constant  pressure 
equal  to  the  each  peak  pressure  for  each  half  period  cycle  and  calculated  the  energy  for 
such  an  alternating  square  wave  signal.  This  was  compared  with  the  energy  calculated 
accurately  for  the  2  /j.  second  recording  using  20  bubbles  as  a  test  population,  we  found 
that  the  accurate  calculation  gave  67  %  of  the  energy  of  the  approximate  square  wave 
calculation.  This  value  is  very  close  the  rms  amplitude  of  each  half  cycle,  which  is  71 
%  of  the  peak  amplitude.  Based  on  this  calculation  we  calculated  the  acoustic  energies 
radiated  by  the  bubbles  by  using  only  the  peak  pressures  of  the  half  cycles  and  taking 
the  67  %  of  the  result  obtained.  For  example,  for  one  of  the  bubbles  we  studied.  plex  , 
the  consecutive  half  cycle  peak  pressures,  were  0.38,  -0.5,  0.23,  -0.19,  0.14,  -0.15  etc. 
pascal  with  an  oscillation  period  of  72  a  sec.    Using  the  formula  : 


Energy  =  (2rr  -j^T )    }  p\ax  At 


where  the  squared  pressures  are  0.14,  0.258,  0.055,  0.04,  0.02  etc.  Pa2  ,  A  t  =  36  n  sec- 

R2 
ond.  R  =  20  cm,  and  2n— — -=  5.66  x  10"  The  sum  of  squared  pressures  is  0.6  Pa2  and 

therefore  the  energy  is  1.06  x  10-12  joule. 

The  average  values  for  the  impact  and  bubble  acoustic  noise  are  shown  in  Fig- 
ures 12  and  13,  with  standard  deviations  and  with  maximium  and  minimum  energies 
obtained  for  three  different  drop  sizes.  The  energy  histograms  for  impacts  and  bubbles 
obtained  during  the  normal  incidence  experiment  are  also  shown  in  Figures  14  and  15. 
The  ratio  of  the  total  average  acoustic  energy  (impact  and  the  bubble)  to  the  calculated 
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kinetic  energy  of  the  drops  just  before  the  impact  is  1.43  x  10~6  for  0.83  mm  diameter 
drops,  6.6  x  10"  for  0.914  mm  and  4.51  x  10-7  for  0.985  mm  diameter  drops.  Where  the 
Kinetic  energies  are  1.58  x  10~6  Joule  for  0.83  mm,  2.74  x  10-6  Joule  for  0.914  mm  and 
3.61  x  10"6  Joule  for  0.985  mm  diameter  drops. 
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Figure  12.  Average  acoustic  energies  radiated  by  the  impacts:  Average,  maximum 
and  minimum  energies  with  standard  deviation  for  normal  incidence 
drops 
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3.     Frequency  Specta 

The  central  frequency  for  the  impact  was  found  by  taking  its  Fourier  Transform. 
A  frequency  spectrum  of  impact  noise  for  a  normally  incidenting  0.83  mm  diameter  drop 
is  shown  in  Figure  16.  Histograms  of  the  impact  peak  frequency  for  0.83  to  0.985  mm 
diameter,  normal  incidence  drops,  falling  with  their  terminal  velocities  are  shown  in 
Figure  17.  Because  the  impact  signal  approximates  a  doublet  (two  opposed  delta  func- 
tions), its  spectrum  is  very  broad.  Its  peak  at  about  15  kHz  agrees  with  the  prediction 
by  Nystuen  (1986). 

The  bubbles  are  similar  to  Type  Al  bubbles  in  breaking  waves  (Medwin  and 
Beaky,  19S9),  which  were  described  as  simply  damped  spherical  bubbles.  The  period  is 
easily  found  by  measuring  the  duration  between  two  consecutive  peaks.  The  frequency 
is  found  simply  by  taking  the  reciprocal  of  the  period.  The  resonance  frequencies  of  the 
bubble  radiation  for  different  sizes  of  drops  are  shown  in  Figure  18. 
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Figure   17.      Peak  frequencies  of  impacts:     Tor  150  normal  incidence  drops  (50  of 
each  size). 
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Figure   18.      Resonance  frequencies  of  bubbles:     For  150  normal  incidence  drops  (50 
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4.     Radiation  Pattern 

To  find  out  the  radiation  pattern  of  the  impact  and  the  bubble  noises  we  first 
used  the  8  hydrophone  array  shown  in  Figure  19,  with  the  hydrophones  placed  every  20 
degrees  on  each  side  of  the  vertical  axis  with  a  highly  sensitive  hydrophone  on  the  center 
for  triggering.  First  it  was  verified  that  the  polar  radiation  pattern  is  nearly  same  as  the 
theoretical  dipole  radiation  pattern,  which  is  the  pgx  cos  6  as  described  in  the  previous 
chapter.  Since  the  pressure  on  both  sides  of  the  axis  was  verified  to  be  symmetric,  the 
array  was  changed  as  shown  in  Figure  20,  with  hydrophones  placed  only  on  one  side  of 
the  axis  with  10  degrees  between  them  to  provide  a  more  accurate  beam  pattern. 

The  polar  radiation  patterns  for  normal  incidence  drops  of  different  sizes  are  as 
shown  in  Figure  21  including  their  theoretical  radiation  patterns.  As  we  can  see  from 
Figure  21,  the  radiation  pattern  is  not  an  exact  dipole  but  is  very  close  to  a  dipole. 
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Figure  21.  Polar  radiation  pattern  of  normal  incidence  drops:  Polar  radiation 
pattern  for  0.3  (i  1  (0.83  mm  diameter),  0.4  n  1  (0.914  mm)  and  0.5  p  1 
(0.985  mm)  drops  at  normal  incidence  and  terminal  velocity  compared 
with  theory. 
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C.     OBLIQUE  INCIDENCE 

The  same  experimental  procedure  used  for  normal  incidence  drops  was  followed  to 
record  and  observe  the  acoustical  properties  of  impacts  and  entrained  bubbles  for  ob- 
lique incidence  drops.  A  fan  was  used  to  impart  a  horizontal  velocity.  To  obtain  the 
angle  of  incidence  and  the  impact  velocity  two  different  methods  were  used.  The  first 
technique  was  to  take  still  photographs  of  drops  in  their  path  by  using  a  stroboscope 
calibrated  to  a  certain  frequency  and  measuring  the  distance  the  drop  travelled  between 
two  consecutive  flashes  of  the  stroboscope.  This  value  was  divided  by  the  time  between 
two  consecutive  flashes  to  find  the  velocity  of  the  drops.  A  5  cm  x  5  cm  string  grid  was 
used  behind  the  path  of  the  drops  to  obtain  maximum  accuracy.  The  angle  of  incidence 
was  measured  by  using  the  same  grid.  The  second  technique  was  to  take  a  motion  pic- 
ture of  the  drops  using  a  video  camera.  Drop  angles  and  velocities  were  then  obtained 
by  using  a  video  recorder  that  has  the  capability  of  slow  motion  (frame  by  frame)  with 
adjustable  speed.  ,\Ye  mostly  used  the  second  technique  which  made  the  observation  of 
the  angle  and  the  impact  velocity  much  easier.  To  measure  the  distances  for  the  second 
technique  again  the  string  grid  was  used. 

1.     Peak  Axial  Pressure 

The  peak  axial  pressure  for  the  bubbles  and  the  impacts  of  different  sizes  of 
drops  were  measured  for  different  incidence  angles.  The  axis  is  again  defined  as  being 
perpendicular  to  the  water  surface.  The  average  pressures  are  shown  in  Figures  22  and 
23.  for  different  drop  sizes  and  different  angles  of  incidence. 
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Figure  22.      Average  peak  axial  pressures  for  impacts:     For  different  drop  sizes  as 
a  function  of  incidence  angle. 
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Figure  23.      Average  peak  axial  pressures  for  bubbles:     For  different  drop  sizes  as 
a  function  of  incidence  angle. 
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2.     Total  Acoustic  Energy 

Total  acoustic  energy  radiated  both  by  the  bubble  and  the  impact  were  calcu- 
lated by  using  the  formulas  described  in  the  previous  chapter  for  a  dipole  perpendicular 
to  the  surface.  The  maximum,  minimum  and  average  acoustic  energies  are  shown  in 
Figures  24  to  26.  The  average  acoustic  energies  radiated  by  the  oblique  impacts  of  dif- 
ferent drop  sizes  are  shown  in  Figure  27.  It  was  observed  that  at  high  incidence  angles 
the  energy  radiated  by  the  impact  increases  as  can  be  seen  from  Figure  27. 

Maximum,  minimum  and  average  energies  radiated  by  oblique  incidence  bub- 
bles for  different  drop  sizes  are  shown  in  Figures  28  to  30.  Average  energies  radiated 
by  bubbles  of  different  drop  sizes  are  as  shown  in  Figure  31.  As  we  can  see  from  this 
figure  average  acoustic  energies  of  bubbles  do  not  change  significantly  with  incidence 
angle.  For  all  these  figures,  around  450  drops  that  created  bubbles  are  used  (around  150 
for  each  size). 
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Figure  24.      Max.  min.  and  avg.  impact  acoustic  energies  for  0.83  nun  drops 
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Figure  25.      Max.  mill,  and  avg.  impact  acoustic  energies  for  0.914  mm  drops 
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Figure  26.      Max.  min.  and  avg.  impact  acoustic  energies  for  0.985  nun  drops 
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Figure  27.      Average  impact  acoustic  energies  of  the  drops  at  oblique  incidence 
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Figure  28.      Radiated  bubble  energy  Cor  0.83  mm  drops:     Maximum,  minimum  and 

average  acoustic  energy  radiated  by  the  bubbles  of  oblique  incidence 
drops,  observed  for  dilfercnt  incidence  angles. 
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Figure  29.      Radiated  bubble  energy  for  0.914  mm  drops:     Maximum,  minimum  ami 

average  acoustic  energy  radiated  by  the  bubbles  of  oblique  incidence 
drops,  observed  for  different  incidence  angles. 
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Figure  30.  Radiated  bubble  energy  for  0.985  mm  drops:  Maximum,  minimum  and 
average  acoustic  energy  radiated  by  the  bubbles  of  oblique  incidence 
drops,  observed  for  different  incidence  angles. 
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Figure  31.      Average  bubble  acoustic  energies  for  oblique  incidence 
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3.     Frequency  Spectra 

The  spectrum  for  the  impact  was  found  by  taking  the  FFT  as  in  the  normal  in- 
cidence case.  The  frequency  spectrum  of  a  typical  impact  at  oblique  incidence  is  shown 
in  Figure  32.  Impacts  of  all  drop  sizes  in  our  range  (0.3  fi  1  to  0.5  \i  1,  or  in  other  terms 
0.83  mm  to  0.985  mm  diameter)  had  broad  spectral  peaks  of  nearly  at  the  same  fre- 
quencies for  all  different  incidence  angles.  The  peak  frequency  shifted  slightly  to  lower 
frequencies  for  oblique  incidence  with  respect  to  the  normal  incidence  case.  The  average 
peak  frequencies  of  the  impacts  is  shown  in  Figure  33. 

Frequencies  of  bubbles  were  found  from  the  period  of  the  bubbles,  following  the 
same  method  as  used  for  normal  incidence. 

The  range  of  resonance  frequencies  of  bubbles  for  0.83  to  0.985  mm  drops  for 
oblique  incidence  are  as  shown  in  Figures  34  to  36.  The  average  frequencies  of  the 
bubble  radiation  for  different  sizes  of  drops  and  for  different  incidence  angles  are  shown 
in  Figure  37. 
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Figure  32.      Frequency  spectrum  of  an  impact  sound:     For  a  0.914  mm  (diameter), 
drop  at  10  degrees  incidence  angle. 
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Figure  33.      Avg.  peak  frequencies  of  impact  sound  for  oblique  inc.  drops 
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Figure  34.      Max.  min.  and  avg.  resonance  frequencies  of  bubbles  for  0.83  mm  drop 
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Figure  35.      Max.  mill,  and  avg.  resonance  frequencies  of  bubbles  for  0.914  mm  drops 
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Fi"ure  36.      Max.  min.  and  avg.  resonance  frequencies  of  bubbles  for  0.985  nun  drops 
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Figure  37.      Avg.  resonance  frequencies  of  bubbles  for  oblique  incidence  drops 
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4.     Percentage  of  Drops  Creating  Bubbles 

The  percentage  of  drops  that  create  bubbles  decreases  sharply  with  increasing 
incidence  angle. This  is  shown  in  Figure  38.  Because  bubble  creation  percentage  de- 
creases with  increasing  angle,  the  ratio  of  the  average  bubble  energy  to  the  impact  en- 
ergy decreases  sharply  as  well,  with  increasing  incidence  angle.  The  ratio  of  the  average 
bubble  energy  to  the  impact  energy  is  as  shown  in  Figure  39. 
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Figure  38.      Percentage  of  drops  that  create  bubbles 
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Figure  39.      The  ratio  of  the  average  huhble  energy  to  the  impact  energy:     The  sharp 

decrease  in  the  ratio  of  the  bubble  enery  to  the  impact  energy  is  mostly 
due  to  the  sharp  decrease  in  the  percentage  of  the  drops  that  create 
bubbles  with  increasing  incidence  angle. 
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5.     Radiation  Pattern 

To  determine  the  radiation  pattern  of  the  impact  and  of  the  bubble  sound  the 
8  hydrophone  array  shown  in  Figure  19  was  used  with  the  hydrophones  placed  even'  20 
degrees  on  each  side  of  the  vertical  axis  with  a  highly  sensitive  hydrophone  on  the  center 
as  in  the  normal  incidence  case.  The  polar  radiation  patterns  of  bubbles  for  different 
sizes  of  drops  and  for  different  incidence  angles  are  as  shown  in  Figures  40  and  41. 
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Figure  40.  Polar  radiation  pattern  of  bubbles  for  10  degrees  incidence:  Radiation 
pattern  of  0.83  mm,  0.914  mm  and  0.985  mm  diameter  drops  for  10 
decrees  incidence 
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Figure  41.  Polar  radiation  pattern  of  bubbles  for  20  degrees  incidence:  Radiation 
pattern  of  0.S3  mm,  0.914  nun  and  0.985  nun  diameter  drops  for  20 
degrees  incidence 
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IV.     CONCLUSIONS 

There  have  been  several  theoretical  studies  and  laboratory  experiments  to  under- 
stand rain  noise  at  the  sea.  Previous  experimental  studies  have  considered  only  normal 
incidence  drops,  often  at  less  than  terminal  velocities. They  have  concentrated  on  either 
bubble  noise  (Crum  and  Pumphrey)  or  impact  noise  (Franz,  Nystuen).  All  the  theore- 
tical studies  have  also  considered  only  normal  incidence  (Longuet-Higgins,  Prosperetti, 
Nystuen).  But  rainfall  consists  of  oblique  incidence  as  well  as  normal  incidence  drops. 
For  this  reason  we  have  studied  both  normal  and  oblique  incidence  drops.  We  have  ob- 
served several  new  and  interesting  results  as  well  as  verified  several  previous  exper- 
imental results. 

We  can  conclude  that  for  both  normal  and  oblique  incidence,  if  a  bubble  is  created, 
the  individual  bubble  has  more  energy  than  the  individual  impact  and  the  total  acoustic 
energy  (bubble  plus  impact)  of  underwater  sound  generated  by  raindrops  is  much  less 
than  the  kinetic  energy  of  the  drops. 

For  normal  incidence,  bubble  noise  increases  with  decreasing  drop  size.  On  the 
contrary,  impact  noise  increases  with  increasing  drop  size.  Both  impacts  and  bubbles 
have  spectral  peaks  around  15  kHz.  The  radiation  pattern  for  bubble  noise  is  very  close 
to  dipole. 

For  oblique  incidence,  the  impact  energy  increases  with  both  increasing  drop  size 
and  increasing  incidence  angle.  When  bubbles  are  formed,  the  energy  per  bubble  does 
not  change  significantly  with  changing  incidence  angle  but  increases  with  decreasing 
drop  size  (as  in  the  normal  incidence  case)  for  all  angles. 

The  percentage  of  drops  creating  bubbles  decreases  sharply  with  increasing  incidence 
angle.  At  20  degrees  of  incidence,  which  corresponds  to  wind  speed  of  1.3  — ,  less  then 
5  %  of  the  drops  create  bubbles  for  the  three  different  drop  sizes  that  we  have  studied. 
Considering  this  sharp  decrease  in  the  percentage  of  drops  creating  bubbles,  we  can 
conclude  that  the  average  acoustic  energy  radiated  by  the  bubbles  decreases  significantly 
with  increasing  incidence  angle.  This  explains  the  reduction  in  the  magnitude  of  the 
spectral  peak  observed  in  natural  rainfall  with  increasing  incidence  angle,  in  other  words 
with  increasing  wind  speed.  The  average  resonance  frequencies  of  bubbles  increase  for 
oblique  incidence.  On  the  other  hand  the  impact  spectral  peak  shifts  slightly  to  lower 
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frequencies.  These  factors  explain  the  broadening  and  shifting  of  the  spectral  peak  ob- 
served in  natural  rainfall  with  increasing  wind  speed. 
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APPENDIX  A.     PREAMPLIFIER  1 

The  block  diagram  of  the  preamplifier  to  obtain  a  fiat  58  dB  amplification  between 
5  kHz  and  45  kHz. 


2.5   N 


Uw 


1    MEG 


0> 


> 


20K 


-VW-j 

200  K 


r 


& 


>  220  K 


>  22K 


2N 


61 


APPENDIX  B.     PREAMPLIFIER  2 

The  block  diagram  of  the  second  type,  simpler  preamplifier  to  obtain  a  flat  55  <JB 
amplification  between  5  and  45  kHz. 
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